Currents in complex polymers: an example of superstatistics for short time series 
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We apply superstatistical techniques to an experimental time series of measured transient currents 
through a thin Aluminium-PMMA- Aluminium film. We show that in good approximation the 
current can be approximated by local Gaussian processes with fluctuating variance. The marginal 
C^l 1 density exhibits 'fat tails' and is well modelled by a superstatistical model. Our techniques can be 

generally be applied to other short time series as well. 
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£ ! I. INTRODUCTION 

a- 

Many time series generated by complex systems in nature are superstatistical, i.e they consist of a superposition of 
^» 1 several dynamics on well-separated time scales. Often, there is locally a simple dynamics (for example, a Gaussian 

process) and the parameters of that simple process fluctuate on a much larger scale. Such varying parameters describe 

i 1 a changing environment of the local system under consideration. Often the relevant measured time series consists 

i-C ] locally of a Gaussian process, with the variance of those Gaussians evolving on a longer time scale. In nonequilibrium 

statistical mechanics, the technique of superstatistics was introduced in [l| and has since then provided a powerful 

Stool to describe a large variety of complex systems for which there is change of environmental conditions [2-9]. A 
superstatistical complex system is mathematically described as a multi-scale superposition of two (or several) statistics, 
one corresponding to local equilibrium statistical mechanics (on a mcsoscopic level modeled by a linear Langevin 
equation leading to locally Gaussian behavior) and the other one corresponding to a slowly varying parameter a of 
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the system. Essential for this approach is the fact that there is sufficient time scale separation, i.e. the local relaxation 



time of the system must be much shorter than the typical time scale on which the parameter a changes. There is 
interesting mathematics associated with superstatistical complex systems. For example, in a recent paper Hanel, 
Thurner and Gell-Mann showed that the superstatistical distribution function f(a) cannot change in an arbitrary way 

■ under macroscopic state changes of the system, and that the underlying transformation group is the Euclidean group 
c"j ' in 1 dimension. 

In most applications in nonequilibrium statistical mechanics the slowly varying control parameter a, underlying the 
superstatistical dynamics, is the local inverse temperature /3 of the system, i.e. a = (3. However, in general the control 
parameter a can also have a different meaning. For example in mathematical finance, a is a vooatility parameter 
7— I ■ describing changing market behavior and turbulences in the market. For a time series, a is usually a local variance 
^ \ parameter that can be directly extracted from suitably chosen time slices. There are numerous interesting applications 

■ of the superstatistics concept to real- world problems, for example to train delay statistics hydrodynamic turbulence 
[lH and cancer survival statistics [Hj]. For further applications, see [T3 - [T9j . 

In this paper, we apply for the first time superstatistical techniques to a complex systems of solid state physics. We 
are interested in the statistical properties of currents flowing through thin films when a voltage is applied. Our mea- 
sured system is a thin aluminium-polymethylmethacrylate- aluminium film (A1-PMMA-A1). A detailed experimental 
investigation of this system was presented in [2(| |2lj j. A particular feature of this system is that when a voltage 
is applied, the corresponding measured currents exhibit transient (nonstationary) behavior with strong fluctuations. 
This can be formally associated with transient chaotic behavior, and in the above papers formally a small positive 
Liapunov exponent was extracted from the data. More recently, Yalcin et al (22l . l23ll made a statistical analysis of the 
fluctuating current data using g-statistics, inspired by the fact that q-statistics [24J is often used for weakly chaotic 
systems which have approximately zero Lyapunov exponent (25l - [27| . 

In this paper we show that the measured current data can be well modelled by a superstatistical dynamics. However, 
in contrast to previous superstatistical time series analysis techniques described in [9J, one has a typical problem for 
these kinds of data: The time series is not long enough to apply the standard superstatistical techniques developed 
in Still, in this paper we show what one can typically do for a short time series. Our method described in the 
following is quite simple and general and can be applied to other small data sets in a similar way. 
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FIG. 1: The current magnitude differences I(t) = i(t + 1) — i(t) for the transient current through thin A1-PMMA-A1 film at 
T = 295K [H 
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FIG. 2: The current magnitude differences I(t) — i(t + 1) — i(t) for the transient current through thin A1-PMMA-A1 film at 
T = 303K 



II. ANALYZING THE EXPERIMENTAL TIME SERIES 



Figs. 1-3 show typical time series of increments I(t) = i(t + 1) — i(t) of the measured current i(t) through the Al- 
PMMA-A1 film at three different temperatures. The average value has been substracted, and all data are rescaled by 
dividing through the standard deviation a determined from the entire time series. We define u{t) = (I(t) — (I(t)))/a. 

One clearly observes regions of strong activity interwoven with periods of much calmer behavior, similar as for 
turbulent flows of share price evolution data. To analyse these data, we divide the time series into a couple of 
windows with qualitatively similar behavior. These windows are shown as vertical lines in Fig. 1-3. 

The choice of the windows is rather arbitrary, they should be large enough to contain enough statistics, but still be 
significantly shorter than the entire time series. 

For each local window, we look at a histogram of the fluctuating currents. As shown in Fig. 4-7, locally the behavior 
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FIG. 3: The current magnitude differences I(t) = i(t + 1) — for the transient current through thin A1-PMMA-A1 film at 
T = 313K 
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FIG. 4: Gaussian fit to the observed local probability density for data points 1-300 of the measured time series of current 
differences I(t) at T = 303K. 



is well approximated by a Gaussian distribution 



P (u\f3) = 




However, the local variance parameter j3 varies strongly from window to window and is thus itself a random variable, 
as expected for superstatistical systems. 



Tab. 1 shows the values of extracted variance parameters /3.; for the various windows and for different temperatures 
T where the experiment was performed. 

We can only extract a few values j3i from our short time series. The general idea of superstatistics is that the 
marginal distribution p(u) is given as a superposition 

p(u) = f f((3)p(u\f3)df3 (2) 
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FIG. 5: Same as Fig. 4 for data points 301-600 of the measured time series 
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FIG. 6: Same as Fig. 4 for data points 601-800 of the measured time series 



for a suitable f(j3) describing the distribution of /3-values. For a short time series, as in our case, one can only make 
guesses of the relevant /(/?)• For this it is useful to look at the histogram of the entire time series, sampling up the 
behavior in all windows. This is shown in Fig. 8. 

Apparently, the histogram data of the entire time series are well-fitted by g-Gaussians of the form 

p(u)~(l + ^(9-l)u 2 )^?. (3) 

The corresponding parameters j3 and q are shown in Tab.l. 

We cannot, as for longer time series extract the relevant distribution f(f3) from a histogram of /^-values, since 
there are only 4 data points. But the fitting result <j3j) can be used to theoretically predict f((3). In our case /(/?) is 
in good approximation a x 2 (or redistribution, meaning that if we were able to look at much longer time series, then 
the /3j would be distributed according to 
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FIG. 7: Same as Fig. 4 for data points 801-996 of the measured time series 
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FIG. 8: Entire PDF of the current magnitude differences I(t) = i(t+l) — i(t) for the transient current through thin A1-PMMA-A1 
at temperature T = 303K [H| 

This follows from doing the integral ([2} using the distribution ((4]), which directly leads to the observed fitting result 
The relation between the parameters n and q is (28j 
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TABLE I: The values of extracted variance parameters /3i, /?2, 03, /34,and /3o for the various slices, for the measured current 
through thin A1-PMMA-A1 film at temperature T = 295, 303 and 313-R". Also shown q and 0. 



FIG. 9: f(/3) of the measured time series through A1-PMMA-A1 film at temperature T = 295, 303 and 313K 



Moreover, the parameter f3 is related to the average value of (3 by [2£ 



(3 := / f(/3)(3d(3 = { ^—^p (6) 



The predicted superstatistical distributions /(/?) are plotted in Fig. 9 for the three different temperatures. We can 
also calculate the average value of Po directly from the time series, namely from the values Pi fitted in the 4 windows, 
i.e. Po — j{Pi + P% + Pz + Pa)- The result is shown in Tab. 1 as well. This sum has huge error terms, because only 
4 windows enter, and the partioning into windows is rather arbitrary. Nevertheless, the obtained average values are 
consistent with the theoretically predicted result © within the statistical error bounds. 

To summarize, for very long time series, one has theoretical methods to decide whether a given time series is 
superstatistical or not [9j. For short time series, as in our case, these methods are not applicable. Still one can do the 
method proposed in the current paper quite generally. First one looks whether a partitioning into a few windows is 
possible such that locally Gaussian behavior is observed. Then one determines the local variance parameters Pi. From 
a histogram of the entire time series one obtains a guess what the relevant distribution /(/3) could be — in our case a 
^-distribution but in general these can be other distributions such as e.g. the lognormal or inverse-;^ 2 distribution 
[2|. After that guess, one can check consistency of derived parameter relations within the statistical error bounds, in 
our case given by relation © but in general given by other relations, depending on /(/?). 
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